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Abstract 

The balance equations of energy-momentum and spin together with Einstein’s field 
equations are investigated by the Liu procedure to find constraints for the consti¬ 
tutive equations in such a way that the Second Law is satisfied. Special cases such 
as spinless systems and curvature insensitive materials are shortly discussed. 

1 Introduction 

The balance equations of continuum physics are formulated for arbitrary materi¬ 
als, that means, they present an underdetermined system of differential equations 
which need for its solution additional equations, the material equations which made 
the system of the balance equations solvable for a special material. There are two 
different procedures: The material equations are presupposecf] and solved together 
with the balances asking later, if the Second Law is satisfied for all positions and 
times. The second method takes the Second Law into consideration before solving 
the balance!! asking for material equations supplementing the balances equations 
so that their common solution is in agreement with the Second Law. 

Here, we deal with the second method which is well known in non-relativistic ma¬ 
terial theory [DEI, but is up to now not strictly developed in general-relativistic 
constitutive theory, except of an attempt [3] • A reason for that delay may be caused 
by the question, how to handle Einstein’s field equations in this framework. Pre¬ 
condition for taking Einstein’s field equations into consideration is that those and 
the balance equations are formulated in a compatible fornE Having done this, the 

*Dedicated to Gerard A. Maugin on the occasion of his 70th birthday 
fCorresponding author: muschik@physik.tu-berlin.de 
lborzeszk@itp.physik.tu-berlin.de, borzeszk@t-online.de 

1 The practitioner knows his material and chooses the correct ansatz. 

2 The theorist looks for a class of materials in agreement with the Second Law. 

3 The energy-momentum tensors are different the in field and balance equations, that needs a modified 
Belinfante-Rosenfeld procedure, here also taking the field equations into consideration [4], 
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question arises, if Einstein’s field equations fit in the well known procedure of gen¬ 
erating the material class being compatible with the Second Lawjh This is the case 
due to the fact that Einstein’s gravitational equations are differential equations of 
second order for the metric. Thus, our approach is completely in accordance with 
Maugin’s opinion [6^1 . Encouraged by his remarks, we derive the Liu equations 
and the dissipation inequality for the combined system of energy-momentum and 
spin balances -that are the Mathisson-Papapetrou equations- and Einstein’s field 
equations. Special cases are considered: spinless systems and curvature insensitive 
materials. 

The paper is organized as follows: After some introductory considerations on tetrads 
and connexions, the balances and the gravitational equations are written down in 
(3+l)-decomposition. State spaces of first and second order are introduced. The 
Liu equations and the dissipation inequality are derived for a first order state space 
material. A discussion finishes the paper. 


2 Tetrads and Connexions 

Measuring values produced by measuring devices need for their description a frame 
of reference which is locally spanned by tetrads {ef}, A,i = 1, ...4, consisting of 
three space-like vectors and one time-like vector numbered bv A These vectors 
represent a standard clock and three standard measuring rodfjj. The field of tetrads 
{ef(x a )} is called an observer field or a system of reference. 

The corresponding contravariant components of the tetrad ef, denoted by e A , are 
defined by 

k _ e/c . . k _ xB /-i \ 

e A e i ~ * -* e A e k ~ °A- l 1 ) 

The tetrad e k A is presupposed to be ortho-normalized, with respect to the metric 
()j k of a pseudo-Riemannian space 

9ik(x a )e i A (x a )e k B (x a ) = r) AB , (r) AB ) := ^ q t J V (2) 

The matrix (tjab) has the shape of the Minkowski metric. Multiplying (j2ib by 
efej? and using ([!]) results in 

ft* = VA B e?e% w f* = if D eZe n D . (3) 

The 16 components { e l A } or {ef}, are restricted by 10 constraints, (J2])i or (]5|). 

Thus for a given metric, 10 components of the system of reference are fixed. 

4 fot the non-relativistic case see [5] 

5 Maugin: “First, from the formal point of view, we may consider Einstein’s theory of general relativity 
as a generalized continuum theory of the second gradient of the space-time metric.... Second, the Maxwell 
stress and other such electromagnetic stress tensors....are none other than purely spatial parts of the space- 
time energy-momentum tensor....Finally, while we can have some doubt with special relativity, it is clear 
that general relativity is a continuum theory from the start. This is more than enough to ponder the 
formulation of continuum thermodynamics in its framework.” 

6 The tetrad representation is well known. For more details we refer to ung. 

7 That is the reason why we introduce tetrads. 
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The tetrads depend on the event, whereas (t]ab) does not. 

The tetrad components e\ and ef, respectively, are matrices transforming 
tensor components at each point of the pseudo-Riemannian space. For a tensor 
of first order, we have the following 

a B := e^dj <—>■ aj = efa B , (4) 

and generally, we obtain for tensors of higher order 

• b... = V- d ..ey B ..., / " 6 ... = T c - D ..e a c e°... . (5) 

We now consider local transformations L B (x a ) of the tetrad components 

at x a : ef = L%(*o) e f, L%(o*)L%(*o) = 5%. (6) 

Inserting (J6]) into (j3j), results in 

a A = L B (*o) a B . (7) 


The properties of the matrix L B are connected to the ortho-normalized tetrads: 
Inserting (J6]) into ([3j) results in 



9ik 


VabLc e ?L 


B % D 
D e k ■ 


( 8 ) 


The Principle of Relativity now asserts that all systems of reference which are 
compatible with the given metric are equivalent to each other, i.e. that e d 
in ([6]) belongs also to an ortho-normalized tetrad as ef in (J3]) does. According 
to (JED, this demand is achieved by the setting 

r] AB Lf(x a )Lf ) (x a ) = r]cD■ (9) 

This equation is the definition for L B (x a ) to be a Lorentz transformation at 
each point ( x a ) . Thus, L B (x a ) belongs to the group of local space-time depen¬ 
dent Lorentz transformations and transforms an orthogonal-normalized tetrad 
into an orthogonal-normalized one at each x a . 

In this paper, we consider the (anholonomic) tetrad representation of Gen¬ 
eral Relativity, instead of the (holonomic) metric one. Accordingly, one has to 
work in a framework, in which the Principle of General Relativity is realized by 
the covariance of the basic laws with respect to local Lorentz transformations, 
instead of general coordinate transformations (what requires the definition of 
the Lorentz connexion)!^]. This is the relativistic version of the Principle of 

8 To avoid misunderstandings, it should be mentioned that the tetrads here introduced as anholonomic 
coordinates differ as well from those introced in tetrad theories of gravitation, where they are considered 
as gravitational potentials, as from those ones introduced as so-called directors for describing the spin of 
the relativistic continua by Maugin and Eringen [F. 
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Material-Frame Indifference (or objectivity)^. 

The covariant derivative dpi is related to the partial one a B ^ m by the con¬ 
nexion as follow^! 

a B\ m = + - n mB a Q, ( 10 ) 


and the connexion Pt® nB is defined by 

n — P A _ri / _|_ n A P Q 

u )in L mi c q ' • 

Here T q mk is the connexion belonging to the metric g ik in (J3])i [IS] 


( 11 ) 


I "mi 2 ^Smq,i T Qiq,m Qmi,q)9 


bq 


( 12 ) 


According to O) i. (TTT1) and (1121b the connexion is a function of the tetrads 
and their first partial derivatives 


q a 


J mQ\ 


= uj“r,(e7 , erf). 


:P 
T >' 




(13) 


3 Balance Equations 

In this section, we start out with the well-known general-relativistic balance 
and field equations of continuum physics. These balance equations are gene¬ 
rated from the special-relativistic balances by applying Einstein’s Principle 
of Equivalence , mathematically realized by the Princple of Minimal Coupling 
which asserts that the partial derivatives in the special-relativistic balances 
have to be replaced by covariant ones with respect to local Lorentz transfor¬ 
mations m- 

3.1 The general case 

Thus, we obtain in a pseudo-Riemannian space the following thermodynamical 
balances: 

Particle number density: 

N A \ a = 0, N a := 1 nu A . (14) 

Here, u A is the material 4-velocity of the particles and c the vacuum velocity 
of light. 

Energy-momentum: 

T AB | fl = f A . (15) 

9 For the special-relativistic theory, being covariant under rigid Lorentz transformations, this was 
already shown by Maugin and Eringen [10]. Later on in a series of papers in JMP EH, Maugin considered 
the space-time covariant description of constitutive equations also in general relativity as the relativistic 
replica of the classical notion of objectivity where covariance means: covariance with respect to general 
coordinate or local Lorentz transformations. For the non-relativistic case see Banana ng. 

10 for more details concerning the tetrad calculus see here and in the sequel [8) 
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Here, T AB is the energy-momentum tensor and f A the external 4-force density. 

Spin: 

(16) 


Sec I A — m BC — —mcB- 


Here, S BC is the spin tensor and m B c the external 4-momentum density. 

Entropy: 

S A lA -<p = a>0 (17) 

with the 4-entropy vector S A , the entropy supply <p and the non-negative en¬ 
tropy production a 0. This inequality represents the Second Law of classical 
held theories, called the dissipation inequality dEum 

Here, the constitutive theory is especially developed in the framework of 
GRT1 12 I. Consequently, we have to take into consideration Einstein’s 

Field Equations: m 


Rab ~ ~VabR — bQab- 

Here, Rab is the symmetric Ricci tensor and R the curvature scalar 


(18) 


Rar R 


c 


lAB lAACBi R Rai 


(19) 


and Radb i s H ie curvature tensor of a pseudo-Riemannian space 


dC 

rt ADB 


f> b P c p m p n pi 

B^mni^b t A t D t B 
-<b 


(r b _ r 6 j_r J h - rt r b i p c p m p n A 

mi ? n L rnn^i ~ L mi L nj L mn 1 ~ ij ) c b C 'A C 'D C 'B‘ 


( 20 ) 


As in the appendix [8] derived, the curvature tensor is linear in the second 
partial derivatives of the tetrads 


r>b _ Tpbuvui ( A\G , mb t A A \ 
±L mni J - J mniG\ c 'j J^u ) vw ' ± mniK^j > 


( 21 ) 


and depends also on their first derivatives and on the tetrads themselves. 
Completely written in tetrad components, we have 


DC _ nuuvw t A A \ Gr I ttU ( A A \ 

ADB — (j ADBG\ t j ) ,vw T Bl ADB \e.j , e ]p ) 


iCuvw 


G 


tC 


( 22 ) 


In accordance with fflSlh the gravitation generating energy-momentum tensor 
Qab has to be symmetric and divergence-free 

©ab = @bu A © AS | S = 0. (23) 

The held equations (TT8|) yield by multiplication with rf :A and by use of (fT9j) o 

R = — k ©b t Rab = b-®ab — ^^ab^Qc- (2 4) 

11 All these quantities are densities. 

12 General Relativity Theory 
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Obvious is, that the gravitation generating energy-momentum tensor ®ab is 
different from that in the energy-momentum balance, because Tab does not 
satisfy (1251) . Therefore the question arises: what is the relation between Qab 
and Ta r' r :i \ An additional question is: are the balance equations (1I5|) and (1I6|) 
compatible with the held equations (TT8ll ? To answer this question, we refer to 
[4], where we proved the following 


■ Proposition: If we presuppose that in GRT energy-momentum and spin sa¬ 
tisfy balance equations - (1151) and (1161) - then the gravitation generating energy- 
momentum tensor in Einstein’s held equations is 


qAB _ rpAB _ 


gCAB gUAU _|_ gAtiU _|_ gt 

and the external force density and the momentum density are 


2“ R’ 

gCAB , qABC , qBAC 


f A = 

m BA = 2 T [BA ]_ 


2 R-cdeS 


CDE 


(25) 

(26) 

(27) 

(28) 


Consequently, the corresponding balance equations of energy-momentum and 
spin are the Mathisson-Papapetrou equations 


rpAB 

T 

gCAB 


1 


_ * rB C' 

A ~ T^^CDE^ 


CDE 


C 


— 2 T^ ab \ 


(29) 

(30) 


The gravitational equations (1241) of General Relativistic Continuum Theory 
(GRCT) become 


R 


AB 


+ yr 4B 'teg = k(t ab - Us CAB + s ABC + s B - ,c )| c ) 
= k(t< ab > ~±(s ABC + s BAC )\ c y 


(31) 


The trace of the gravitation generating energy-momentum tensor is according 
to (13TT) and (124|) 

(32) 


f\D _ rpD nDC 

~ 1 D \q- 


Consequently, the Ricci tensor in (l3Tj) becomes 

1/^/1 nr* A r 1 ' 1 


R 


AB 


= K (T^-^S ABC + S BAC ) lc -^ V AB 


rpD _ q_ 

1 D J D I c 


~iDC 


(33) 


which is determined by the symmetric part of the energy-momentum tensor 
and derivatives of the spin tensor. The held equations (l3Tj) become 

(. S ABC + S BAC + 7 1 AB SR c ) lr , + -R ab = 2T (as) - r) AB TE- (34) 

l G K 
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For the Special Relativity Theory, this question is answered in m- 
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According to (l34jh all derivatives appearing in the field equations are on the 
lhs. This form of the field equations is essential for the exploitation of the 
dissipation inequality (1171) by the Liu procedure which is discussed in sect El 


3.2 (3+l)-Decomposition 

For the exploitation of the dissipation inequality, we need (l29]h (130|) and 
(1341 in (3+l)-decomposition. Starting out with the particle flux, the energy- 
momentum tensor, the spin tensor and the 4-entropy [231 241 [25] , we obtain 


N a = \nu J 
c z 

4 

(35) 

T AB = (i e «8 + ip B )^ + i tI V + U B , 

(36) 

S CAB = ( -L- 4 ® + i- u ' A ~. B ')v. c + t-u lA ~ B V + s CAB , 

\C Z C 4 / c z 

(37) 

S A = \su A + s A . 
c z 

(38) 

In the sequel, we split the fields of the (3+1 ^decomposition into two classes. 
Using the projector h AB onto the sub-space perpendicular to the 4-velocity u A 

h AB ;= AB _ 1 n A u B = bBA^ 

c z 

(39) 

the first class is defined by 



n 

= N a u a , 

(40) 

e 

= TabU A u b , £ := e/c 2 , 

(41) 

sab 

= S BF liAh B uc, 

(42) 

Sn 

= 2S% F u c u E h F A . 

(43) 

Here n is the particle density , e the energy density , sab the spin density and 

finally the spin density vector. 


The second class of fields of the (3+l)-decomposition is defined by [23, 25] 

t AB 

= h AE T EF h FB , 

(44) 

P A 

= h AE T FE u F , 

(45) 

d A 

= h AE T EF u F , 

(46) 

sab 

= S EF liAh B hQ, 

(47) 

77 A 

^B 

= 2S EF hQU E h B , 

(48) 

S 

= S A u a , s A := S B h A . 

(49) 

Here t AB is the stress tensor , p A the momentum flux density , q A 

the energy 

flux density , s^ B the couple stress, E B the spin stress, s the entropy density 


and s A the entropy flux. 


7 




We now write down the balance equations in (3+l)-decomposition. According 
to (125|) and (1361) . we obtain the energy-momentum balance 


i(pjX + u B qff) + tf B + ^(p B uf A + uf^) + i(«V) |, 


— dB qCDE 

— 2 rLcDEj 


(50) 


The spin balance results from (1301) and (1371) 


-u 


l a £ b]c + S CAB + — u [A Z B]C + 


AB 


+ ^U 


[A^B]\ U C 


i\o 


= 2 T [ab \ (51) 


I o I C C 2 | C L \ C* 

and from (1771) and (1381) follows the dissipation inequality 

^ s \a uA + S\a + ^ su \a - V- ( 52 ) 

The field equations (134)) by use of (1371) result in a long-winded expression 


1 


-U 


[■ B^C]A + g ABC + _ yf B ~p\ A + 


C 


\c 


2“|W 


+ -^ B ~ C 1 
L V C l C 4 


U 


i\c 


+ 


+ \u^ ]B + s BAC + \u [ ^ B + [ (\s AC + \u^)u B 


c 


C 


:; 2 \c 


\C 


+ 


+ rr{j 2 u {D ^ + S g° |c+ ±u IDlc E^+ [(X s g + 


: C]D DC 


\C 


+ 


-R ab = 2 T {ab) - r) AB T%. (53) 


The system of diffential equations (I5H to (1551) is valid for arbitrary materials. 
For solving it, we need constitutive equations and the Second Law (1521) has to 
be taken into account. 


4 State Spaces 

4.1 Basic and constitutive fields 

The balances (fill) . (1T51) and (TTfil) and the field equations (TT8|) represent equa¬ 
tions to calculate the 58 fields N A , T AB , Sf c , S A and ef by taking the con¬ 
straint (jSJ) into consideratiorJ 44 !. Additionally, we have a further constraint by 
the dissipation inequality (1171) . Balances and field equations are valid for arbi¬ 
trary materials, e.g. non-Newtonian fluids, spin-fluids, solids, liquid crystals. 
To solve the system of differential equations combined of the balances and the 
held equations, we have to introduce constitutive equations. Therefore, we 


14 4+16+24-|-4-|-10=58. The 10 gtk are replaced by 10 bilinear combinations of the 16 ef according to 











have to split these 58 fields into those which represent the independent vari¬ 
ables, called the basic fieldlfi\, and those which are material dependent ones, 
named the constitutive fields [26]. 

Like in non-relativistic continuum physics, we choose as basic holds the first 
class (HUD to (HUD and additional the held of tetrads which according to ©1 
are tensors of hrst order in the tetrad indices 

z = (n,u A ,e,s A B,Z A ,ef) =: ( w,ef ). (54) 

These basic variables do not contain any derivative. Taking the normalization 
of the 4-velocity 

u a ua = c 2 —» u a ua\b = 0 (55) 

and (J3]) into account, the number of the basic helds in (154|) is 21@ 

According to the splitting into basic and constitutive helds, we now have to 
introduce 58 - 21 = 37 constitutive fiehk0. These are 

M=(t AB jP A ,q A ,s c AB ,E A ,S A ). (56) 

Basic and constitutive helds are tensors of different order, that means, they 
are covariantly defined. 

4.2 Constitutive mappings 

For describing constitutive properties, we have to introduce material map- 
pina i^l which connect the constitutive helds (HUD . At", with the basic helds 
(151) . z, thus characterizing the material. Obvious is, that the basic helds 
as the independent variables do not describe materials sufficiently in non¬ 
equilibrium. Consequently, we have to extend 2 by additional variables £ 
consisting of derivatives of the basic helcl^HI 

C (^\mi Z\m\m • • ■) • (57) 

The order of the additional derivatives which are enclosed in £ depends on 
the considered material. Consequently, the set of constitutive properties (J56D 
depends on (z, £). which is called the covariantly dehned state spac^l and 
which spans the domain of the material mapping A4 

M = M(z, C). (58) 

15 often named the wanted fields 

16 1+3+1+3+3+10=21 

17 9+3+3+9+9+4=37 

18 also called constitutive mappings m 

19 An other possibility of extending the basic variables is to take instead of derivatives constitutive fields 
into the extension, e.g. t AB and q A , see: Extended Thermodynamics [281129] , 

20 or constitutive space 
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Its range M is spanned by the constitutive properties (1561) . The entire equa¬ 
tion (158|) is kown as the constitutive or material equation 

AT Ad (z, 2 | m , Zj m j n ,...). (59) 


If we use partial derivatives instead of covariant ones, the connexion has to be 
added to the additional variables 


c 




n B n B 
1 114 mAi s 


and the constitutive equation becomes [27] 


M 


Ad(.z, z ^ ra , z^ mn . 


n B ct B \ 


(60) 


(61) 


According to (1T31) . the connexion itself depends on the tetrads and their first 
partial derivatives which are included in the (z !m , Z }Tn ^ n ,...). Consequently, 
4 can be cancelled in the set of th state space variables 


AT — Ad (z, z jm , z mtn ,...). 


(62) 


In this form of the constitutive equation, the state space does not consist 
of tensors in contrast to (159|) . but it contains all non-tensorial quantities so 
that by combination of them, the tensorial constitutive equation (T59|) can be 
rediscovered. Here, we prefer the non-tensorial representation of the constitu¬ 
tive equation, because curvature and Ricci tensor are functions of state space 
variables which according to (1221) are partial derivatives. 


4.3 First and second derivatives state spaces 

Taking (TTT]) o into account, the additional variables (1571) of first derivative state 
spaces are 


C ' 7 


(^ ? m) 


e A 

11 5 r 



(63) 


Consequently, the state space of first partial derivatives becomes 

z 1 = (z,z 9 m) = (w,w 9m ,ef,ef 9 m) = 

(TIE'LL , Sab, i SAB^rm i — A^rm 9 m)* (64) 

The state space Z 1 is a local object defined at x a , spanned by independent 
variables -the wanted basic fields and their first derivatives- generating locally 
the constitutive holds (1561) by the material mapping (158|) . 

The additional variables of second order are 



(- z ,m; 2 , mn ) — ( w ,mi ' w ,mni e i imi e i ,mn) — 

(C ) ^.mni ^A^mni SAB^mm , mn)i 


(65) 
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resulting in the state space of second order 


11 { k _ _i A _ _i A 

I 71, W , e, j ^A^mi ^^mi ^AB^mi 'A^mi ,m> 

^jmn.) ^A^mni £jinni ®AB ^mm ^A^mni ,mn^ • (66) 


Having introduced the concepts of state space and material mapping, we are 
able to consider the system of differential equations (130]) to Q53j) on the chosen 
state spaces. 


4.4 Balances and field equations on the state space 


We especially consider a state space of first order (jMj) and write down the 
system m to (153]) in which expressions like H c B | 4 appear. These expressions 

have according to ((4|) and (1T01) the form 

□ cb m = 0 CB |,4 = (□“ +nf„0° B +(if Q a c “y A . (67) 

Taking (1131) into account, we obtain 

B cs | A = 0“4 + { HS S }(*, (68) 

{□S B }(z,2| m ) := (nf o 0« fl +!Jf Q 0 c «)4. (69) 

Because of (1131) and the chosen state space of first order, the bracket symbol 
{ } does not depend on higher derivatives. The partial derivative has 

the following meaning, distinguishing between basic and additional variables 
according to (153]) 




CB 


ll 


(dB CB k da CB k da CB Q d^ CB 0 

V^F w ” + a^ w « + ~^f et » + 54r y* « 


(70) 


By use of (1701) . (1051) and (1331) . (1001) results in the energy-momentum balance 
on the state space 
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The lhs of the energy-momentum balance (1711) is linear in the higher deriva¬ 
tives w k , m j and whereas no higher derivatives appear on the rhs which 

depends on the state space variables of the first order state space (1M1) . 

Starting out with (152]) . the dissipation inequality becomes with (j68j) 

+ sA ,i e A + M> > ~ ^ SU \ A ' (72) 

By taking d70]i into account, we obtain the the dissipation inequality on the 
state space 

1 ds A ds A \ i k f 1 ds , ds A 

~ 2 o^ u + d^r AW ^ + h^ u + ^ 

> --T-—« A + — 

vc 2 dw k dw k ) A ,J vc 2 de k de k 

+<P ~ { si } - J 2 SU \ A - ( 73 ) 

The pretty long-winded energy-momentum balance (1711) and the dissipation 
inequality (1731) can be written in a symbolic form 

energy-momentum: u A Bmj w k ^ m j + 12 Ag kmj e k ^ m j = 1 C B . (74) 

dissipation inequality: w k + 2 Bq^ e k > V. (75) 

The coefficients 11 A, 12 A, 1 C, '£>, 2 B and V can be directly read off from (ITT]) 
and (1731) . 


7 Q 

e A e k ,3 


+ 


A cQ 

c A c k ^ 


The spin balance (1711) and the field equations (1771) can be treated in the same 
manner as the energy-momentum balance and the dissipation inequality. Their 
symbolic shape is as follows 

spin: 21 A£ Bmi w\ mj + 22 Af tmi e'i, ml = 2 C AB , (76) 

field equations: 21 A ABmj w k tm j + 22 Aq Bkm2 e® — 2 C AB . (77) 

The definitions of 21 A 22 A. 21 A 22 A , 2 C AB and 2 C AB can be obtained from 
(EDi and (1731) by the same treatment which was performed with the energy- 
momentum balance. We will not perform this procedure here, because our 
main aim is to demonstrate that an exploitation of the dissipation inequality 
is also possible in General Relativity Theory. 


In the sequel we need the balances of energy-momentum and spin, the gravita¬ 
tional field equations and the dissipation inequality - (1771) to (1771) - in a common 
matrix formulation 


/ llA B mj 12 jBkmj \ 
21 a ABmj 22 A ABkmj 
31^ABm.j 32 j^ABkmj 


1 gmj 2 jgkmj 
k Q 


W 


•>rn] 


'k 


W 


•>rn] 


Q 

-k irnj 



> V. 


(78) 

(79) 
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The dissipation inequality (T79l) enforces that the A, B , C and V are not ar¬ 
bitrary. The exploitation of the dissipation inequality, the Liu procedure, is 
considered in the next section. 


5 Exploitation of the Dissipation Inequality 

The usual way to introduce the Second Law -the dissipation inequality (TT71) - 
into the exploitation of the balance equations -here (TT4jh (l29]h (1301) and (134)) - is 
to solve this system of coupled differential equations and afterwards to check, 
if the dissipation inequality 11171) is satisfied for all events x a . Apart from 
the complexity of this procedure, it is not possible, because for solving the 
balances, we need constitutive equations determining the constitutive fields 
(l56j) . But here, we are not interested in solving these balances, but we are 
looking for general restrictions to the constitutive equations by the dissipation 
inequality. 

5.1 The Liu procedure 

On account of the unknown constitutive equations, we cannot solve the bal¬ 
ance equations, but we can ask for the conditions which must be satisfied by 
the constitutive equations, so that the dissipation inequality is locally valid 
for all events. The Liu procedure is such a tool for exploiting the dissipation 
inequality with regard to the balance and field equations resulting in restric¬ 
tions for the constitutive equations OG0 The procedure removes the higher 
derivatives which are outside of the state space, and it generates constraints 
for the derivatives of the constitutive properties (156|) in the state space. Here, 
we present the Liu procedure in a symbolical formalism. 

The balances of energy-momentum (1291) and of spin (15U|) . the field equations 
(134|) and the dissipation inequality (TT6j) have the following shape 

K AC( A M) = m?(A M ) > ®. (80) 

Here, A AI is the symbol for the state space variables (z, £) in (15B|) . 

Starting out with (l80]b . we obtain according to (168|) 

= ^ffe j c + {M$ c } = (81) 

We now split the state space variables 

A m = (A^,Af) (82) 

into those ones A+ whose derivatives A + j, the higher derivatives, are out 
of the state space. Then, we split (18T| into that part which is linear in the 

21 An other procedure to take the Second Law into consideration is that of Coleman and Noll US- 


13 




higher derivatives, and into the remaining one which forms a part of the rhs 
of the following equation 


jdM AC (A M ) aM 
c dA f + 


, d^A™} M 

c dam - 


{®g c }+0 A 


m A . (83) 


Analogously, the dissipation inequality (l80j) o results in 


j d ffl c (A M ) 
dA 


A 


M 


• > — p J 

iJ — e C 


dm c (A M ) 
dA F 




= : ©■ 


(84) 


Now we formulate 


■ Liu’s Theorem [1, 2]: There are functions of the state space variables A^(A a/ ) 
generating the Liu equations 


A A ^ C 


m 


aAC 


= e J 


; da 


c 


dA f dA+ 

and the reduced dissipation inequality 

A A a A > ©, 

removing the higher derivatives A + j. 


A. 


dl 


aAB 




B 


dA 


M 


<9Af 


(85) 


( 86 ) 


Liu equations and the reduced dissipation inequality represent constraints for 
the constitutive quantities KI‘ 4C , 0^ and 0 which we derive in the next 
section. 


5.2 Liu equations and reduced dissipation inequality 

According to (j531) and (jM|) . the lhs of the matrix equations (1751) and (1791) are 
linear in the higher derivatives. Consequently, we can apply the Liu procedure 
which according to (18511 ■> and (jSSj) results in the following matrix equations 


( X A b 2 Aab 3 A ab ) 


/ 11 A Brn j 
21 j^ABmj 


31 A ABmj 32 a A 

\ 1/1 k 


12 j^Bkmj \ 


( 3 A b 2 A ab 3 Aab 



> V. (88) 


Thus, we obtain the following Liu equations 

l B™ j = l A B u A Bmj + 2 A ab 21 At Bmj + 3 Aab 31 A£ Bmj , (89) 

2 B k Q mj = X A B 1 2 A Bkmj + 2 A ab 22 A ABkmj + 3 \ AB 32 A ABkmj , (90) 


and the reduced dissipation inequality 

1 A b 1 C B + 2 A ab 2 C ab + 3 A ab 3 C ab > V. (91) 
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According to (174)) to (177)) . lm A belongs to the energy-momentum balance, 29 A 
to the spin balance and 39 A to the gravitational field equations. The same is 
valid for the 9 C. This division allows a classification into special cases which 
some of them are discussed in sectEJ 

5.3 Constraints of state space variables 

The particle number balance (1141) 

n \ a = 0 = (nu A ) | A = ny A u A + rm^ A (92) 

and also (1551) represent according to (jMj) a constraint for the state space vari¬ 
ables. How to handle these constraints with respect to the Liu procedure, is 
expressed in the following 

■ Proposition [5T, :52If55|: Constraints of the state space variables do not influ¬ 
ence the Liu procedure, that means, the results achieved by the Liu procedure 
are independent of introducing the constraints before or after performing it. ■ 


Consequently, the particle number balance and the normalization of the 4- 
velocity do not generate additional constitutive restrictions. 

6 Special Cases 

6.1 Ignoring the field equations 

The Mathisson-Papapetrou equations (1551) and (1551) are compatible with Ein¬ 
steins field equations (1511) . Ignoring the field equations means that the influ¬ 
ence of the material on the geometry is not taken into account: the curvature 
in (12 9 1) is given ad-hoc, but compatible with the field equations. Decoupling 
the field equation means, we have to ignore in (157)) and (1SS1) all quantities 
having the index combination 3, @. The Liu equations (1891) and (155)) become 

'h? = 1 As 11 A^ mj + 2 A A B 2l Af m \ (93) 

2 Bq J = 1 An 12 Ajf m V l AR 22 A A Q Bkm \ (94) 

and the reduced dissipation inequality is 

1 An 1 C b + 2 A ab 2 C ab > V . (95) 

The constraints on the entropy - (1931) and (1941) - and on the entropy produc¬ 
tion (155|) are changed with respect to the case, that the field equations are 
taken into consideration - (1551) to (1511) . Ignoring the field equations results in 
altered material properties, if the Second Law is taken into account by the Liu 
procedure. 
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6.2 Vanishing spin 

If the spin tensor is set to zero, we obtain from (127)) . (1251) and (157)) 

T [ab] = 0, T^ B) = 0, ^ R ab = 2 T {ab) - ri AB T%. (96) 

The spin balance is satisfied by (l96lb and (]96]) o. Consequently, we ignore all 
quantities with the index combination 2, @. The energy-momentum balance is 
GU with the restriction that only the symmetric part (AB — » ( AB )) comes 
into play and that S CDE is set to zero. Thus, 11 A, 12 A and l C can be read off 
from (171)) . 

The gravitational held equations (]96)L become with (122)) 

2 -G D A%i Q 4, mi = -Hg DB + 2T {AB) -r, AB Tg. (97) 

rt 

To take into consideration the Second Law, the following matrix elements 
vanish for state spaces of first order 


21 .A, 22 ^l, 2 C, 31 A are zero, 

(98) 

and we obtain 


o 

32 a kmj ^ /-yDkmj 

^ ABQ ~ ADBQ 

r\j 

(99) 

3 Cab = —H b db + 2Tp4s) — VabT b . 

(100) 

Consequently, the Liu equations (1891) and (1901) become 


l B^ = 'A B ll A* mi , 

(101) 

2 jgkrnj _ 1 12 j^Bkmj ^^ B 32^ABkmj 

(102) 

and the reduced dissipation inequality is 


1 k B l C B + 3 A ab 3 C ab > V. 

(103) 


This example of the Mathisson-Papapetrou equations demonstrates again that 
the dissipation inequality and its exploitation have to be supplemented, thus 
taking the Second Law into account. Even in the easy case of spinless material, 
the Second Law cannot be ignored independently of taking the gravitational 
held equations in the course of the Liu procedure into account or not. 

6.3 Curvature insensitive material 

A material is called curvature insensitive, if its material properties do not de¬ 
pend on In that case, the material does not contribute to the curvature, 
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because the last term vanishes in (170]) and with it the influence of the mate¬ 
rial upon the curvature. Of course, the type of curvature sensitive material 
appear, always if the curvature comes into play. 

The curvature does not explicitly appear in the spin balance (130]) and in the 
dissipation inequality (1171) or (107]) . Consequently, the matrix elements 22 A 
and 2 B vanish identically for curvature insensitive material according to (1751) 
and (1701) . According to ([90]) . we obtain 

o = !A fl ^Ag kmj + 3 Aab 32 AQ Bkmj , (104) 

an expression which also for curvature insensitive materials depends on the 
curvature because the curvature appears explicitly in the energy-momentum 
balance (170]) and in the field equations (13TT) . 

The entropy production (101]) may be influenced by the curvature: if the A 
depend on the curvature, the entropy production depends on it, too. 

7 Discussion 

The balance equations of continuum physics are formulated without taking a 
special material into consideration, that means, one needs constitutive equa¬ 
tions to obtain by them a system of differential equations which is not under¬ 
determined. Constitutive equations are not arbitrary: they must be modeled 
in such a way, that the solution of the resulting system of differential equations 
satisfies the so-called material axioms. Among other axioms, the Second Law 
may be the most prominent one. In continuum physics, it runs as follows: 
Constitutive equations must have the property that the entropy production 
is never negative for all times and positions. A tool for finding such suitable 
constitutive equations is the Liu procedure which is here performed in the 
framework of General Relativity Theory. 

The characteristic of general-relativistic continuum physics is, that beyond the 
balance equations, Einstein’s gravitational field equations have to be taken 
into ccount. In accordance with Maugin’s view described in footnote 0, we 
formally consider Einstein’s gravitational equations on the same footing as the 
balance equations of continuum thermodynamics, rendering to apply the Liu 
procedure for the system of equations consisting of the energy-momentum and 
spin balances ~ (fT5l) and (fl6l) ~ together with Einstein’s equations (fT8l) . 

Compatibility of the balances with the field equations enforces, that the sources 
of the energy-momentum balance and of the spin balance are not arbitrary, 
but result in the Mathisson-Papapetrou equations, independently of including 
the field equations into the Liu procedure or not. As expected, the constraints 
for the constitutive equations are different, if the gravitational field equations 
are included into the Liu procedure or are ignored. Two other special cases 
are shortly discussed: the spin-less system in General Relativity and curvature 
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insensitive materials which are the only materials in Special Relativity, but a 
special case in GRT. 

In the present paper, the Liu procedure is performed for material mappings 
defined on a first order state space. This results in the Liu conditions — ff89|) 
and (M|- that must be satisfied in order to guarantee the compatibility of 
the combined system of differential equations - (TT5jh (HRl) and with the 

Second Law of thermodynamics (TT71) . 

Due to the complicated structure of the Liu equations in GRT, it is extremely 
difficult to arrive at general statements with regard to the compatibility of 
the general-relativistic continuum mechanics with the Second Law. Starting 
out with a constitutive equation by a substantiate guess, the Liu conditions 
- fl89]i to fl9Tf - differentiate between gravitational fields as solutions of Ein¬ 
stein’s equations: there could exist solutions for which the constraints on the 
constitutive equations can be satisfied and such ones for which this cannot be 
achieved. Regarding the compatibility of the solutions of Einstein’s equations 
with continuum thermodynamics as an essential physical requirement, the Liu 
conditions do not qualify all solutions of the gravitational equations as to be 
physically relevant. In this sense, the Liu conditions act as a censor for the 
constitutive equations with respect to the Second Law. 

Because the covariant curvature tensor (1221) cannot be formulated covariantly 
as a function of the state space variables (see: sect JS]) , consequently the Liu 
procedure cannot be performed in a covariant manner. That is no disadvan¬ 
tage bearing into mind, that special solutions of the gravitational equations 
cannot generally be characterized covariantly, either. 


8 Appendix 

Curvature tensor in tetrad representation 

Starting out with (J3])i, we obtain from (1T2|) that the connexion with respect 
to the metric is a function of the tetrads and their Erst partial derivatives 


V b (p F p F ) = 


= 2 WT+ ef , m ) + effef - ef , ff ) + ef (ef, m - e*, q )}e b D e%105) 
We obtain from (fT21) by partial differentiation 

r mi,n — 7) {dmq^in + 9iq,mn ~~ 9mi,qn)9 9 + ^ mi9pq9 9 ,nj (10®) 


resulting in 


pO _ pO 

1 mi^n 1 mnj* 


~~ 2 ,mn Qmi^qn Qnq^mi + 9mn,qi)9 ? T 9pq^X mi9 ^rnn9 9 ,i)- (107) 
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Taking ©i and 


Qmn^qi VAb(& 



<£„) (108) 


into account, we obtain that the curvature tensor depends linearly on the 
second derivatives of the tetrads. Consequently, it has the form which is given 
in (Ell) . Some easy, but long-winded calculations result in 



(109) 


( 110 ) 


The final result is generated by inserting (l3lb and (1105)) into (II 101) . 
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